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1 Introduction 



Naked singularities [T] have been objects of great interest in the theory of general 
relativity for some time now. With a proof of the cosmic censorship hypothesis 
still not in place, attempts have been made to understand possible realistic physi- 
cal scenarios involving these objects. For example, the work of [2] models galactic 
centers with naked singularities, and studies them in the context of gravitational 
lensing (see also [3]). On the other hand, efforts have been made to confirm 
that these are stable objects, and there are hints from numerical work that this 
might be the case |1] (more recent attempts to settle the question of stability 
of naked singularites were made in [5]). Recently, the Janis- Newman- Winicour 
(JNW) naked singularity [6], [7] has also been shown to be stable under scalar 
field perturbations [8]. 

It is important to gain knowledge about naked singularities, and to contrast 
them with their black hole cousins. Several works in this direction appear in the 
literature, but to the best of our knowledge, there has so far been no attempt 
to calculate radiation spectra in naked singularity backgrounds. In this paper, 
we study the scalar power radiation spectrum in the background of the JNW 
naked singularity. Power spectra in black hole backgrounds have been studied in 
several works (see, e.g [9], [10], [H]) using perturbation theory. Here, following [9] 
(see [12] for a more recent treatment in the context of AdS black holes), we 
consider a scalar field perturbation [8J of the JNW naked singularity, with a 
massive particle coupled to this field, and calculate the scalar power radiated. 
Mostly in the literature, such work has been done in connection with synchrotron 
radiation, i.e where most of the radiated power lies in a range of frequencies that is 
large compared to the orbital frequency of the particle. In fact, the work of [9j was 
initiated with the hope that such radiation would be observable, although later 
work [T3] suggested limitations on the astrophysical implications of gravitational 
synchrotron radiationjl] Here, we present a comprehensive analysis of the scalar 
power spectrum emitted by a particle in a circular orbit, coupled to a scalar field, 
in the JNW background. We do this both for stable as well as unstable orbits, by 
employing a parabolic WKB approximation. We find that for a certain range of 
parameters, the radiation spectrum is different from what is obtained with black 
hole backgrounds. For a different range of parameters, there are two distinct 

^In a slightly different context, radiation from particles in orbit around black holes has been 
investigated in a series of papers [14] . 
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types of spectra, one of which mimics a black hole. 

This paper is organized as follows. In section 2, we review the physics of 
the JNW metric and set up the basic equations for calculating the scalar power. 
Section 3 discusses the methodology and our main results, and we conclude in 
section 4 with some discussions. 



2 The JNW Metric and the Scalar Equations 



In this section, we start by reviewing the basic features of the JNW metric and 
circular geodesies therein, following the work of [15]. We then proceed to set 
up the equations for scalar radiation in the background of this metric, which will 
enable us to calculate the power spectrum in the next section. Using the notation 
of [7], we start with the JNW metric [6] : 



Here, v varies between and 1. It can be verified that this space-time has no 
horizon. Further, it can be shown [17] that it satisfies the weak energy conditions 
and that the singularity is globally naked. The JNW space-time is sourced by a 
scalar field. 



where q is its magnitude. The ADM mass of the space-time M is related to 



eq.([I]) is given by z/ = 2M/h. The Schwarzschild limit of the JNW singularity is 
obtained by taking z/ = 1, i.e g = 0. We will specialize to the case where we have 
a particle moving in a circular orbit in the JNW geometry, setting 9 = 7r/2, using 
the spherical symmetry of the system. The geodesic equations appear in [15] and 
we briefly review them. For time-like geodesies, the equations following from the 
metric are. 





(2) 



the parameters b and q hj b = 2^ (f ^ M^, and the parameter v appearing in 




(3) 



h 



(4) 




f2 + W{r) = 



(5) 
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where, 



W{r) 



1 



1 



-7 



(6) 



The dots denote differentiation with respect to an affine parameter, and h and 7 
are conserved quantities, corresponding to the angular momentum and the total 
energy. The angular frequency of the particle is 

(7) 



h 



coo = — = — 1-- 
dt r'^'j \ r 



For a circular geodesic, f = 0, and hence from eqs. ([5]) and ([6]) 



7' = W{ro) 



To 



(8) 



where tq denotes the radius of the orbit. After imposing -7^ 



^l^^^jj = 0, we obtain 
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^2ro-bu-b) 
(2ro - 2biy - b) 



h = To 



bu 1-^ 



2rn -2bv-b 



(9) 



From eq.Q, we see that the radius of the photon sphere [2], where the energy 
and the angular momentum diverges, is at 

b 



ps 



2(1 + 2^) 



(10) 



Depending on whether the photon sphere exists (i.e is above r = 6) or not, 
the singularity is called a weakly naked or strongly naked singularity [2J. From 
eq. lfTOl) . we see that for u > 1/2, the JNW solution becomes a weakly naked 
singularity, while it remains strongly naked for < z/ < 1/2. The quantity Vps 
will be important for our discussion on the power spectrum. From eq.Q, we get 
an useful expression for the angular frequency of the particle : 



bu 



^ N (2u-l)/2 

To J \2ro-bu 



1/2 



It is also necessary for our purpose to derive some restrictions relating to the 
radius of the orbit, tq. Such restrictions for stable circular orbits have been 
worked out in [TJ] and we briefly state their results. Using the notation 

b 



r± = -{l + ?>u± V5z/2 _ 1 j 



(12) 
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for < < for a stable circular orbit, we require 

6 < ro < oo (13) 

When < < |) there are two different regions where stable circular orbits 
exist: 

6 < ro < r_ and r+ < ro < oo (14) 
Finally, for ^ < z/ < 1, stable circular gedesics exist for 

r_|_ < ro < oo (15) 

For future reference, we record the expression for m° = 

V roj \2ro-2bu-bJ ^ > 

Now, in the JNW space-time, we introduce a scalar field $, that does not 
interact with the field if) (sourcing the fixed JNW metric), and treat $ as a 
perturbation over the JNW background. That the JNW metric is perturbatively 
stable under such a scalar field has been demonstrated in [8] . A massive particle is 
coupled to the field and we study scalar emission from this particle. To describe 
the physics, we write down, following [9] , the action describing the interaction of 
the particle with the field $, in the JNW background, 

5 = y d^x^ $,^<l>'^ -moj dT{l + g,$) v^=i^ (17) 

Here, z^{t) denotes the world line coordinates of the particle, its scalar charge 
and mo its mass. The relevant part of the energy momentum tensor is, 

T^"' = ^ (^^'^'^'-^ - ^g"''^^^^^'^^ + mo(-^)^/^ j dr^l + qs^)5\x - z)z^'z'' 

(18) 

The equation of motion for $ is 

^99"'^A=^^qsmA — ] 6\x - z{t)) (19) 



dx^^ \ dx^ J \ dr 

The delta function in the preceding two equations can be expanded as [9] : 

6^{x - z{t)) = 6{r - ro)6{e - 7r/2)5(0 - ooot) (20) 
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Using the spherical symmetry of the JNW background, we can write 



oo oo 



m=—oo l=\m\ 



-imwQt 



(21) 



and consequently, we obtain the radial equation for uimir) as. 



drl 
where, 



rn^ujQ I 1 



Cim = -47r(M°ro)-ig,moFr(^/2,0) 
and the Regge- Wheeler coordinate r* is given by; 

— = !-- = /(r) 
dr^ r 



uim = CimS{n-ro*) 
(22) 

(23) 
(24) 



We will focus on the solution of eq.( l22|) . and use it in conjunction with the power 
radiated, given by [9] 



mod / dnr'r: = — I dnr^^^ 

Att I ot or 



(25) 



Having elaborated on the basic setup, we now proceed to calculate the scalar 
power spectrum in the JNW background. 



3 Scalar Power Radiation Spectra in the JNW 
Background 

First, we find the solution for the homogeneous part of the radial equation of 
eq. fl22|) . Then, the full solution can be found by a Green's function technique. 
Following [8], the asymptotic behaviour of the homogeneous part of the radial 
equation is, 

uim ~ Ar^ + B; r^, -oo (26) 

and 

uim ~ Ce~"* + De-~''* ; r, ^ oo (27) 

Here, A, B, C and D are constants. Let ul and ur be the left moving and 
right moving solutions of the radial equation. We can then write the boundary 
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conditions 

f (A;oo)-'/'(e-*^'-* +5e*'^*); n ^ +00 
Ml ~ < (28) 

I An + B; r* — —00 

Ur ^ < (29) 
I Cn; n — )■ — cx) 

where 5 and T are analogues of reflection and transmission coeflicients, and k^o 
is a normalization factor. The Wronskian is computed as 

W{+oo) = -2iT, W{-oo) = -BC (30) 

The Green's function is, 

G(r,,roJ = { ' (31) 

From eqs.f l28p and (!29|) . the main difference of the asymptotic solutions here, as 
compared to the black hole case can be seen. Namely, near r = b, the solution 
is not a plane wave. A similar situation occurs in the case of AdS black holes, 
near spatial inflnity p^, and the wave function is set to zero there. That the 
situation demands care in the choice of initial conditions as r — 6 is of interest 
by itself, and deserves further study. We will conflne ourselves to the calculation 
of the power radiated towards inflnity. Using standard properties of the Green's 
function of eq.( l3T]) in eqs.llTTl) and ( l25l) . the expression can be shown to be [9] 

^ 00 00 

^=^5^ Vmu;o|Gm|VL(ro*)r (32) 

m=-oo l=\m\ 

Next, we attempt to solve the radial equation, and hence flnd an analytic solution 
for M^, by a WKB approximation. Speciflcally, we use the parabolic WKB method, 
described in [9]. We will very briefly recapitulate the method here, and point out 
a few issues regarding its apphcation in our case. Let us suppose we have an 
equation of the form, □ 

+ - V] y{x) = (33) 



^In the present discussion, we use the variable x to keep the discussion general, following [5]. 
We keep in mind that x will be taken as r* in our case. 
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Further, we assume that the potential can be represented near its maxima xq by a 
parabola, V{xo + x) ~ V{xo) + ^V" {xo)x'^, with V" (xq) < 0. Using the parabohc 
approximation for the potential in eq. fl33p . the general expression for the wave 
function in the region near the maximum of the potential barrier is given by [9] 

/ 1 r(-z/) /I 1 1 1 \ , , 

where 

w=\^V"{xo)\^^; e = ^^"^"j^"^ ; = (1 - = -1(1 + ^e)(35) 

Dy{rix) is the parabolic cylinder function. We choose the potential maxima in the 
parabolic WKB method to be near xq = 0. In defining the tortoise coordinate, 
we have the freedom to choose the integration constant such that the potential 
maxima does indeed occur near tq* = 0. Then, in eq. fl32l) . we can use 



where a standard expression for -0,^(0) has been used. Of course, one has to 
be careful in applying the results of [9] to the present problem. Recall that the 
parabolic WKB approximation relies on the fact that given a potential which can 
be adequately represented by a parabola near its maxima, there exists classical 
turning points {V = E\ and that there is an overlap region, where the parabolic 
solution and the standard WKB one are supposed to hold simultaneously. In such 
a situation, the asymptotic solution of eq.( 15^ is matched with a WKB solution. 
In the present case, this poses a problem. Although classical turning points exist, 
from eq. fl28|) we see that asymptotically near the horizon, does not have the 
form of a plane wave, i.e is not of a WKB form, whereas the parabolic cylindrical 
functions asymptote to plane wave solutions in this region. As alluded to before, 
this is a peculiarity in dealing with a naked singularity, which we have to live 
with. In the region — )■ oo, however, such a problem does not arise, and can 
be obtained by matching the parabolic solution with a plane wave [9]. We will 
proceed, with this in mind, i.e our is obtained only from the properties of the 
solution at spatial infinity, and we are unable to comment on its behavior near 
r = h. 



•^An alternative method would be to obtain a numerical solution for with the bound- 
ary conditions specified in eq. ((28|) . as was done for the AdS Schwarzschild black hole in |12) . 
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Since we now have the expressions for the wave function, using eq. (l36|) and 
eq. (!23|) as inputs in eq.( l32|) . we can calculate the power radiated. Of course, we 
also need the expressions for uq and u", given by eqs. lfTTl) and (fT6l) . respectively. 
In our case, we take the Schro dinger- like equation to be of the form of eq. (|33ll . 
with E = 0. That is, the potential is taken (from eq. fl22|) ) as 

\ r J \r-^ / V ^/ 

with uo calculated from eq. (fTT]) . The analogy with quantum mechanical scatter- 
ing theory is not clear, but we are able to solve the equation nonetheless. We 
note here that there is an immediate caveat in using the potential of eq. (137|) in a 
WKB scheme, namely, the peak of the potential, and the value of Uq is dependent 
on the azimuthal quantum number, m. We thus have to resort to a reasonable 
approximation, so that the parabolic WKB method can be trusted. This can be 
done only if the peak of the potential does not shift appreciably, on varying m 
and /. We will justify this in sequel. 

Before we start our analysis on the power spectrum, let us first state what we 
expect. From a gravitational lensing perspective, it can be shown (see, e.g. [3] 
and references therein) that for z/ > |, the JNW naked singularity presents 
qualitatively same features as a black hole. Physically, we expect our analysis of 
the power to retain the same features, i.e for this range of u, the radiation should 
be synchrotron in nature, (i.e, peak at large m), and moreover, the radiated 
power should sharply peak in a narrow band of frequencies - for black holes, it 
is well known [9j that at the peak of the spectrum, more than 99 percent of the 
power is radiated in the / = m mode. However, as we will far as the 

radiated power is concerned, there are interesting deviations from this picture, 
and that only for z/ > | does the power spectrum have the same qualitative 
features as a Schwarzschild background. A word about the nature of orbits is 
in order. We will, in sequel, consider both stable and unstable circular orbits in 
the JNW geometry. The work of [9] specifically focuses on unstable i.e highly 
relativistic orbits in the Schwarzschild geometry, where the radiation was shown 
to be synchrotron. That such orbits might not be of relevance in astrophysical 
scenarios was pointed out in [13j. In what follows, we will not be concerned with 

Unfortunately, however, we found that in our case, such numerical solutions generated with 
standard Mathematica or Python numerical routines could not be trusted beyond moderately 
large values of I and m. 
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Figure 1: V{r) of eg. (1371) vs r with u = 
0.03. The sohd blue, dashed red and dot- 
dashed green hues are for m = 1,1 = 1, 
m = 5,/ = 5 and m = 100,/ = 100 
respectively. 
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Figure 2: Power radiated in the mth 
mode with z/ = 0.03. The dashed red 
line is for the / = m mode, while the 
solid blue line is the power, where I has 
been summed from m to m + 2. 



this latter fact, and our purpose would be to distinguish between black hole and 
naked singularity backgrounds, as far as radiation properties are concerned. 

To begin with, we consider the case < u < We first choose z/ = 0.03, 
so that we are far from the Schwarzschild limit. Recall that for this value of u, 
all circular orbits are stable, from eq. (fT3|) . We will henceforth choose 6 = 2, so 
that the corresponding Schwarzschild mass is set to unity, in appropriate units. 
We will also set to unity the mass mo and the scalar charge of eq. llTTl) . First, 
in eq. fl371) . we substitute Uq from eq. fllll) . Then, in order for the parabolic WKB 
approximation to remain valid, we require the particle orbit to be close to the 
peak of the potential, i.e where ^1^=^^ = 0. While this condition is satisfied for 
two positive values of tq, one of these is always at a position ro < 2, and hence 
discarded. In the present case, we find that for (m, /) = (1,1), V{r) peaks at 
tq = 2.88, while for large values of m = /, the peak asymptotes to tq ~ 2.97. The 
shift in the peak with increase of m, / is thus small, and we choose our circular 

''in the treatment of it is not possible to analytically deal with stable orbits in the 
Schwarzschild background, since the analogue of the potential of eq.(|37|) in that case peaks 
close to r — 3M, the location of the photon sphere, whereas stable orbits in the Schwarzschild 
geometry occur beyond r = 6M , M being the Schwarzschild mass. In our case, it is possible to 
understand radiation from stable orbits. Expectedly, as we will see, this radiation has a large 
angular spread. 
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orbit at tq = 2.92, and apply the parabolic WKB approximation. In fig.(IT]), 
we have shown the potential of eq.( l37|) . for u = 0.03. In this figure, the solid 
blue, dashed red and dot-dashed green lines correspond to {m,l) = (1, 1), (5,5) 
and (100, 100) respectively, where we have appropriately scaled V{r), so that the 
three cases can be displayed on a single graph. It can be seen that the peak of 
the potential shifts towards r ~ 2.97, as we increase the values of m = /. In 
fig.(|2]), we show our result for the power radiated, as a function of the azimuthal 
quantum number m. The dashed red line correspond to the power radiated in 
the m = I mode, whereas in the solid blue line, we have summed over / from m 
to m + 2. It can be seen that the radiation is dominated by the small m modes. 

Here and later in the paper, we focus on the modes close to / = m. This is 
justified by taking a concrete example. In the present case, i.e for u = 0.03, given 
a fixed value of m, the peak of the potential of eq.(|3711 asymptotes to ~ 2.99 when 
I is made very large, so that our approximation of tq = 2.92 can still be applied 
with reasonable confidence. Specifically, for I = m = 2, the potential peaks at 
ro = 2.93, while for (/, m) = (2, 3) and (2, 4), the potential peaks at vq = 2.96 and 
2.98 respectively. Now, for I = m = 2, the numerical value of the radiated power 
is Pi=m=2 ~ 9.4 X 10"^, while P;=4,m=2 ~ 1-7 X 10"^ (Pi=4,m=3 = from properties 
of spherical harmonics). Thus, the radiated power decreases away from the / = m 
mode, and to a good approximation, it is enough to consider modes near / = m, 
to understand the nature of the spectra. We will keep this in mind in our future 
discussions. 

We now come to the region < u < 0.5. Here, we will take u = 0.47 as 
our representative value. We choose vq = 2.7 in our analysis below, and it can 
be checked that the parabolic WKB approximation can be well justified for this 
choice, for small values of m and /. This will be enough for us, as the radiated 
power is negligible for large values of m and /. Note that from eq.f ll2p . this means 
that our particle is in the unstable region, as stable orbits exist for tq < 2.087 
and ro > 2.733, and we have kept the particle very close to the boundary of 
the unstable region. Here, we find that for the I = m modes, the radiated 
power maximizes at small values (~ 4) of m. The main features of the radiation 
spectrum is similar to the previous case, < z/ < i.e it is dominated by the 
small m modes, and we will not elaborate on this further. 

^In hindsight, we note here that when ly = 0.5^, an interesting possibihty occurs, the 
radiation becomes synchrotron, and can resemble those with black hole backgrounds. This will 
be clear from our discussion below. 
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Figure 3: Power radiated as a function Figure 4: Power radiated as a function 

of m for V = 0.6, with tq = 0.6. The sohd of m for u = 0.6, with tq = 2.201. For 

blue curve is the case where / has been the solid blue curve, I has been summed 

summed from m to m + 2. The dashed from m to m + 2. The dashed red curve 

red curve is the I = m mode. is the / = m mode. 



While the cases considered till now correspond to strongly naked singularities, 
more interesting cases occur for 0.5 < z/ < 1, for which a photon sphere exists. 
As an illustration, let us take z/ = 0.6. Here, the radius of the photon sphere 
is = 2.2, and from eq.( !T2|) . = 3.69 and r_ = 1.91. Thus, stable circular 
orbits will exist only for tq > 3.69. In this case, as we show below, there are 
two distinct peaks of the potential, both in the unstable region, one of which is 
very close to the photon sphere. Specifically, for z/ = 0.6, the potential of eq.f l37|) 
extremizes, for m = / = 2, at tq = 2.70 and at tq = 1.8, the latter being discarded 
(we remind the reader that we are interested in the modes near m = I). As we 
increase to m = 400 (=/), there are two maxima, at tq = 2.41 and at tq = 2.19. 
Further increasing the value of m, we find that one of the maxima asymptote to 
rg = 2.4 and the other merges with the radius of the photon sphere. We can thus 
consider unstable orbits very close to rps. 

Now, for small m and Z, we choose our orbit at tq = 2.6, where we expect 
our parabolic WKB approximation to hold good. In fig. ([3]), we have shown the 
spectrum as a function of m, and see that the peak of the spectrum is at m ~ 7. 
In this figure, the solid blue curve denotes the case for which / has been summed 
from m to m + 2. The dashed red curve is the / = m mode. It can be checked 
that taking tq ~ 2.4, i.e close to one of the asymptotic values of the potential 
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Figure 5: Power radiated as a function of pigure 6: Power radiated as a function of 

e in the (small) m = / mode, for z/ = 0.03 q for large m, for v = 0.6 (dotted black) 

(dotted black), 0.47 (dashed red) and 0.6 and 0.9 (sohd red), and for u = 0.6, with 

(solid blue). small m shown in dashed blue. 

maximum, does not change the result appreciably. Fig.(j3]) shows the spectrum 
for To = 2.201, where the same color coding as in fig. ([3]) has been used. The peak 
of the spectrum has changed drastically from the previous case. This, in fact, 
closely resembles a black hole spectrum, where it is known that unstable orbits 
exhibit synchrotron radiation [9]. 

As we increase the value of z/, we find an additional feature in the power 
spectrum. Namely, at u = 2/3, the two peaks of the potential for large values 
of m, as alluded to above, merge. Specifically, for u = 2/3, for I = m = 2, 
the potential of eg. ( 1371) has maxima at tq = 2.71 and 1.89, with the second one 
being unphysical. As we increase the value of m (= /), the two maxima merge at 
ro = 7/3. For u > 2/3, we find that there is a single maximum of the potential in 
the region r > 2, which asymptotes to r^s as the values of m and / are increased. 
While the radiated power maximizes for small to moderate values of m for orbits 
away from rps in these cases, § it shows synchrotron behavior when the particle 
radius is close to rps. As expected, in the limiting case — )■ 1, the potential 
always maximizes close to rps, for all values of m and /. 

We also comment upon a further distinct feature of the power spectrum. To 
illustrate the point, we take u = 0.6, and choose ro = 2.201, close to the radius 
of the photon sphere. Here, the peak of the spectrum occurs at m = 3669. We 

^We always stay within the vahdity of the parabohc WKB approximation. 
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find that the ratio of the / = m and the I = m + 2 modes is ~ 4. In contrast, 
for u = 0.9, choosing tq = 2.801 (rps = 2.8) the peak power is for m = 1448, and 
the same ratio yields a value ~ 24. Hence, as we approach the black hole limit 
z/ — )■ 1, most of the power is radiated in the / = m mode, consistent with the 
results of |9], but away from it, although the radiation might be synchrotron in 
nature, a few modes away from / = m cannot be entirely neglected. 

Finally, we have also studied the angular dependence of the power spectrum, 
at its peak, for various values of u. The results are summarized in figs. ([5]) and ([6]). 
In fig. ([5]), we consider the m mode for which the power is maximized, for small 
m, and / = m. We see that the radiation is spread over a wide angular region, 
for z/ = 0.03 (dotted black), 0.47 (dashed red) and 0.6 (solid blue). In contrast, 
the dotted black and solid red lines in fig. ([6]) show the angular dependence of the 
maximally radiating mode for large m (i.e when the particle is close to the photon 
sphere), for u = 0.6 and 0.9 respectively. We see a relatively small angular spread 
of the spectrum. To make the contrast clearer, we have plotted, in fig. ([6]), the 
small m case of z/ = 0.6, shown in dotted blue. 

4 Discussion and Summary 

In this paper, using a WKB approximation scheme following |9j, we have made 
a comprehensive investigation of scalar radiation in the background of the JNW 
naked singularity. Our main conclusion here is that the nature of the power spec- 
trum reveals a rich structure, for different ranges of the parameter u appearing 
in the JNW metric. For ^ < the spectrum shows a peak at small values of 
the azimuthal quantum number. The situation is unchanged for < < |- 
However, from our discussion, it follows that in this region, close to z^ = 0.5~, 
the nature of the spectrum may become synchrotron. The most interesting case 
occurs for i < z/ < 1, where we saw that upto ^ = ^, there are two maxima of the 
potential of eg. (1371). Thus, there might be two distinct types of radiation spectra, 
one of them closely resembling that from a black hole, when the particle is close to 
the photon sphere. Beyond this value of u, the potential has a single maximum, 
which stays close to the photon sphere as z/ — t- 1. We have also seen that away 
from the z/ — )■ 1 limit, when the radius of the particle orbit is close to that of the 
photon sphere, although the qualitative behavior of the spectrum may be similar 
to a black hole background, there is an important difference. Namely that a few 
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I modes away from I — m, contributes significantly to the power, in contrast to 
the black hole case, where more than 99 percent of the power is radiated in the 
I — m mode. 

The nature of naked singularities is far from being fully understood, and we 
believe that our results complement existing ones in the literature, that study 
the similarities and differences between these objects and black holes. It will 
be interesting to study the nature of electromagnetic or gravitational radiation 
spectra in the background of the JNW singularity. We leave this for a future 
publication. 
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